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Abstract 

The minor phases of powder metallurgy disk superalloy LSHR 
were studied. Samples were consistently heat treated at three 
different temperatures for long times to approach equilibrium. 
Additional heat treatments were also performed for shorter times, 
to assess minor phase kinetics in non-equilibrium conditions. 
Minor phases including MC carbides, M23C6 carbides, M 3 B 2 
borides, and sigma were identified. Their average sizes and total 
area fractions were determined. CALPHAD thermodynamics 
databases and PrecipiCalc ®, a computational precipitation 
modeling tool, were employed with Ni-base thermodynamics and 
diffusion databases to model and simulate the phase 
microstructural evolution observed in the experiments with an 
objective to identify the model limitations and the directions of 
model enhancement. 

Introduction 

Structure, composition, fractions, and size distributions of carbide 
and boride phases can significantly influence the processing and 
mechanical properties of nickel-base superalloys [1-6], especially 
at elevated temperatures. MC carbides are stable up to very high 
temperatures near superalloy solidification, and therefore can help 
constrain grain growth during solution heat treatments above the 
gamma prime (y ') solvus [7], The elements Ta, Ti, Nb, and W 
often combine with carbon to form this carbide. M13C6 carbides 
are stable at lower temperatures up to near 900°C, and can 
improve the resistance to grain boundary sliding, rupture and 
dwell crack growth at service temperatures of 600-815°C. Cr, 
Mo, and W can react with carbon to form this carbide. Boron in 
solid solution within the y matrix also improves the resistance to 
grain boundary sliding, rupture, and dwell crack growth [5]. The 
M 3 B 2 boride phase has not been widely shown to specifically 
impact mechanical properties, but is often observed in disk 
superalloys containing boron. This phase can have a solidus 
temperature only slightly higher than solution heat treatment 
temperatures [6], and therefore can limit solution heat treatment 
upper tolerance temperatures. Mo, Cr, and W often combine with 
B to form this phase. 

Harmful topological close packed (TCP) phases can also form in 
disk superalloys [1,2]. These phases such as a, p, Laves, and P 
have low inherent ductility due to limited slip systems and often 
have needle or lath morphologies, and can sometimes reduce 
tensile and rupture ductilities at application temperatures [8, 9]. 
Several models have been developed over the years to predict the 
formation of such phases [10,11], 

Recent advances in multi-component thermodynamic modeling, 
based on CALculation of PHAse Diagrams (CALPHAD) [12,13] 


framework, have resulted in software packages such as Thermo- 
Calc [14], JMatPro [151 and PANDAT [161. When used in 
conjunction with Ni-base databases, these software tools can be 
used for the prediction of relative percentages and compositions 
of phases of Ni-base superalloys in equilibrium at various 
temperatures. Further development of CALPHAD-based 
precipitation kinetics modeling has extended the capability to 
predict the transient precipitation microstructure with software 
tools such as PrecipiCalc [17,18], JMatPro [15], PanPrecipitation 
[16] and TC-PRISMA [19]. Model calibration and validation of 
CALPHAD databases and models with experimentally measured 
quantities [18,20] is an integral part of integrated computational 
materials engineering (ICME) discipline [21], 

In the first half of this article, measured minor phase identities, 
fractions (contents) and sizes in disk superalloy LSHR are 
reported after known heat treatments. In the second half, the 
results of CALPHAD-based equilibrium predictions made using 
several multicomponent thermodynamics databases and 
PrecipiCalc precipitation simulation are presented. The primary 
objective of this work was to determine the predominant minor 
phases in the disk superalloy LSHR after varied heat treatments, 
and to assess existing model capabilities and limitations for 
predicting these minor phases. 

Experimental Characterization of Minor Phases 

Disk superalloy LSHR having a composition in weight percent of 
3.46 Al. 0.028B, 0.029C, 20.7Co, 12.52Cr, 0.07Fe, 2.73Mo, 
1.45Nb, 0.01880, 0.0014N, 0.03Si, 1.6Ta, 3.50Ti, 4.33W, 
0.049Zr, bal. Ni and trace impurities was produced through 
powder processing. It was atomized in argon, hot compacted, 
extruded and then isothermally forged into a flat, uniform disk 
5.08 cm thick and 30.48 cm in diameter. 

Phase extraction samples were removed by electro-discharge 
machining pins of about 10 g weight and 10 mm diameter in the 
circumferential direction from the rim of the disk. These samples 
were solution heat treated at 1199 °C for 1 h and furnace cooled at 
an average rate of 2 1 °C/min. They were then aged at 843 °C for 
1,000 h and water quenched. Electrochemical extractions of 
minor phase were performed using a 10 ml HC1-1 g tartaric acid- 
90 ml methanol electrolyte. Triplicate extractions were performed 
in each case. Extracted phases and lattice parameters were 
identified using x-ray diffraction (XRD) [22], 

Smaller pins of 3 mm diameter and 16 mm length were also 
machined in the circumferential direction near the rim of the disk, 
for performing various assessments. Samples were extracted from 
the extrusion before forging, and from the forging, for 
examination without further material processing. Heat treatment 
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conditions of the remaining samples are summarized in Table 1. 
Single samples were only solution heat treated at 1199 °C for 1 h, 
then water quenched at an average cooling rate of 5,000 °C/min. 
or furnace cooled at a controlled average cooling rate of 21 
°C/min. Additional samples were solution heat treated at 1199 °C 
for 1 h, water quenched at an average cooling rate of about 5,000 
°C/min to minimize phase nucleation and growth during the 
cooling process, and subsequently given aging heat treatments. 
They were aged at temperatures of 760, 843. and 927 °C for times 
of 10, 100, and 1,000 h and again water quenched, yielding nine 
different aged cases. 


Table 1. Heat Treatment Conditions. 


Sample ID 

1 199 C / 
lh / 

Solution 

Quench 

Aging 

Temp 

-C 

Aging 
Time - 
h 

Aging 

Quench 

LAF-Forged 





LFF 

Furnace 




LFW 

Water 




LFW760-10 

Water 

760 

10 

Water 

LFW760-100 

Water 

760 

100 

Water 

LFW760-lk 

Water 

760 

1000 

Water 

LFW843-10 

Water 

843 

10 

Water 

LFW843-100 

Water 

843 

100 

Water 

LFW843-lk 

Water 

843 

1000 

Water 

LFW927-10 

Water 

927 

10 

Water 

LFW927-99 

Water 

927 

99 

Water 

LFW927-lk 

Water 

927 

1000 

Water 


Selected heat treated and quenched samples were sectioned and 
conventionally metallographically prepared. The sections were 
etched using a waterless Kallings etch. Linear intercept grain size 
was measured from optical images according to ASTM El 12. 
SEM evaluations were used to select samples for determination of 
minor phase compositions using TEM and microprobe. Thin 
disks of 3mm diameter and -0.5 mm thick were sliced, 
mechanically ground, and polished from selected heat-treated 
samples. These disks were then electrochemically thinned to 
electron transparency using a solution of 10 % Perchloric and 90 
% Methanol cooled to -30 °C. A FEI CM200 Transmission 
Electron Microscope (TEM) operating at 200 kV was used for 
imaging and to identify crystal structure of predominant minor 
phases using selected area electron diffraction patterns (SAEDP) 
and energy dispersive X-ray spectroscopy (EDS). Associated 
bright field images were used to estimate the typical size range 
and morphology of each phase. Sections having larger size minor 
phase particles were usually selected for estimation of minor 
phase compositions using a JEOL JXA-8200 scanning electron 
microprobe. At least 5 particles of each minor phase were 
analyzed. 

Sizes and area fractions of various second phases were quantified 
via SEM using metallographically prepared sections, in order to 
sample sufficient areas. Samples were imaged in a Hitachi S- 
4700 Field Emission SEM at 20 kV using a working distance of 
12 mm. Secondary electron and backscattered electron images 
were obtained at 2,000 x magnification, in order to include 4-8 
grains and their associated boundaries within each image. This 
enabled statistically valid area fraction estimates for each phase in 
each image, even for phases such as M 2 jC 6 observed only at grain 


boundaries. Five adjacent images were acquired for each sample, 
with no bias of imaging towards a certain region of particles or 
grain boundaries. An ED AX CDU LEAP™ detector in the SEM 
was used for semi- quantitative chemical analysis of each particle 
visible in these images. Compositions determined from EDS 
spectra were then compared to results previously obtained in thin 
foil TEM evaluations. While the EDS spectra acquired using 
TEM and SEM were not identical, due to different interrogated 
volumes and instrumental conditions, relative major elemental 
compositions for a given phase were consistent for both cases. 
Based on theses analyses, the phase identity of each particle was 
recorded in the SEM images, and later used for size 
determinations of various phases. The combination of these two 
techniques enabled discrimination and measurement of minor 
phase contents and particle sizes over a much larger imaged area 
than what was feasible using TEM. 

Particle sizes were measured using Sigma Scan Pro 5.0 image 
analysis software. The difference in brightness between minor 
and y-y ' phases in secondary electron SEM images enabled 
brightness-based local thresholding. Each minor phase was 
measured separately, and also classified by location, either within 
grains or at grain boundaries. Area, equivalent radius, major axis, 
minor axis, and perimeter were determined for each particle. The 
equivalent radius calculation assumed a spherical shape and used 
the formula: 

Equivalent Radius= Area / n 

Phase area fraction was determined by combining the areas of all 
particles of an identity and class within an image. Size-frequency 
histograms were generated for each phase with more than 10 
particles acquired. 

Results and Discussion 

Heat Treatments 

A typical SEM image of a sample extracted from the forging is 
shown in Figure 1 . Linear intercept grain size of this sample was 
5.5 pm. The large dark phase particles are primary y with mean 
radius of about 1.1 pm. These particles acted to constrain grain 
growth during extrusion and forging. Samples subsequently 
solution heat treated at 1199 °C for 1 h then furnace cooled and 
water quenched are shown in Figures 2 and 3. The primary y ' 
particles have dissolved during this solution treatment. These 
supersolvus solution heat treatments therefore allowed grain 
growth to a larger linear intercept grain size near 25 pm. Only 
MC carbides and M^B^ borides were identified in these samples. 



/ 


M3B2 ^ 

K703 20.0kV 1 1 8mm x2.00k SE(M) 6/2/2009 11,11 'gp^ 1 


Figure 1 . Image of LAF, as-forged sample. 


894 




Figure 2. Image of LFF, forged sample solution heat treated 
1 199 °C / 1 h / furnace cooled. 



Figure 3. Image of LFW, forged sample solution heat treated 
1 199 °C / 1 h / water quenched. 


Images of samples aged 1,000 h at 760, 843, and 927 °C are 
compared in Figures 4, 5, and 6, respectively. MC and M 23 C6 
carbides, M 3 B 2 borides, and sigma ( 0 ) phases were identified in 
these conditions [ 16] . These phases were also identified in 
samples aged shorter times of 10 and 100 h at 760, 843, and 927 
°C. 

Subsequent measurements of minor phase sizes over multiple 
SEM images indicated the equivalent radius of each minor phase 
could be considered normally distributed, so that size distributions 
could be represented using mean and standard deviation for each 
phase. 




1856 20.(*V 12 0mm x200k SE(M) 6/^2009 


Figure 4. Image of LFW760-lk, forged sample solution heat 
treated 1 199 °C / 1 h / water quenched, aged 760 °C / 1,000 h. 



Figure 5. Image of LFW843-lk, forged sample solution heat 
treated 1 199 °C / 1 h / water quenched, aged 843 °C / 1,000 h. 



Figure 6. Image of LFW927-lk, forged sample solution heat 
treated 1 199 °C / 1 h / water quenched, aged 927 °C / 1,000 h. 

MC Carbides 

MC carbides had an FCC Fm3m crystal structure, with a 
lattice parameter of a = 4.3820 angstroms (A) measured with 
XRD. This was comparable to that reported for several other disk 
superalloys [23]. but slightly higher than that reported for Rene 
88DT, an alloy somewhat similar in composition to LSHR [24], 
MC carbide average equivalent radius did not exceed 0.27 pm 
(Figures 1-3), making quantitative microprobe measurements of 
composition difficult. However, EDS measurements of these 
carbides in TEM thin foils indicated this phase had an 
approximate composition of 26Ti-12Nb-12Ta-50C in atomic 
percent, Table 3. Quantitative compositions are not often reported 
for minor phases such as carbides and borides existing at such 
small sizes, due to difficulties of in-situ measurements using 
scanning electron microscopy. However, surveys indicate many 
other disk superalloys also have MC carbides which contain Ti, 
Ta, and Nb [1,2, 23, 24], The slightly higher lattice parameter 
compared to Rene 88DT could be related to the presence of Ta in 
MC carbides for LSHR. Overall MC carbide size, morphology, 
and area fraction appeared comparable in most samples, 
indicating high overall phase stability for the various conditions 
examined. While most appeared spherical or slightly elliptical in 
the SEM images presented. TEM examinations indicated the MC 
carbides usually had short cylindrical or blocky shapes, each made 
up of a single crystal. 

Yet, some changes in MC carbides were observed when 
comparing those within grains versus at grain boundaries. About 
35 % of the carbides were observed at grain boundaries after 
forging. Figure 1 . However, no more than 23 % were observed at 
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grain boundaries for samples subsequently water quenched and 
aged, Figures 2-3 and Table 2. Here they more often appeared 
within the grains, sometimes loosely distributed in clusters. 
Additional systematic evaluations at lower magnifications would 
be necessary to fully quantify the degree of clustering in their 
spatial distribution. The MC carbides at grain boundaries were 
often slightly larger than those within grains. MC carbide content 
at grain boundaries decreased with increasing aging time at 843 
°C, which will be further considered below. 

M 23 C 6 Carbides 

M 2 jC 6 carbides had an FCC Fm3m crystal structure, with lattice 
parameter estimated using TEM of 10.6 A. This is very close to 
that reported for Rene 88DT [24], The average phase 
composition was measured by microprobe on a sample aged 843 
°C/ 1,000 h, having particles sometimes as large as 0.5 pm radius, 
Fig. 4. Microprobe quantitative measurements indicated this 
phase had an approximate composition of 33.0Cr-20.6Co-14.3Ni- 
5.4Mo-2.5W-24.2C in atomic percent. Table 2. Surveys indicate 
other disk superalloys also have M33C6 carbides which are 
principally high in Cr and Mo [1,2, 23, 241 . The high Co and Ni 
content measured here implies these elements can substitute for 
Cr in LSHR. The M 2 3C 6 carbides populated this alloy as particles 
often elongated along the grain boundaries. However, larger 
particles sometimes resided at triple points. TEM evaluations 
indicated the elongated particles were often made up of multiple 
grains. These carbides were not found away from grain 
boundaries. While some of these particles were highly elongated, 
equivalent radius was calculated to enable consistent comparisons 
among all particles and minor phases in this study. 

Mean equivalent radius and area fraction of M 21 C6 carbides are 
shown versus sample condition in Table 2. These carbides were 
not found in any samples that were not given aging heat 
treatments. This indicates that significant time is required at the 
aging temperatures for carbide formation to occur, exclusively at 
grain boundaries. Results show consistent trends in growth and 
dissolution for the different times and temperatures investigated. 
At both 760 °C and 843 °C, the average equivalent radius and area 
fraction both increase with time. This trend was enhanced at 843 
°C, where the longest aging time of 1,000 h produced highest 
phase area fraction of 0.0042 and equivalent mean radius of near 
0.29 pm for this phase. Coarsening occurred at this temperature 
beyond 100 h, where number density decreased while particle 
radius increased. Recall from Table 2 that MC carbide content at 
grain boundaries decreased with increasing time at 843 °C. No 
MC carbides partially transformed in-situ to M23C6 carbides were 
directly observed for these samples, as previously reported for 
several other nickel-base superalloys [1,2], Due to the difficulties 
in isolating such particles in SEM of metallographic sections and 
TEM of thin foils, this transformation could still be possible here. 
Based on the much wider spacing of prior MC carbides compared 
to that of M23C6, it appears that nucleation and growth of new 
particles accounted for a large majority of M23C6 phase particles 
observed in aged specimens. No M 2 3C 6 carbides were observed 
within grains in all aged conditions, and many MC carbides were 
still present there. At 927 °C, M33C6 was found after 10 h aging. 


but not for longer times, indicating it is not stable at this 
temperature. Thus as the aging time increased and conditions 
approached equilibrium, the size and area fraction were driven to 
zero. 

M , B , Borides 

M,B 2 borides had a tetragonal P4/mbm crystal structure, with 
lattice parameters estimated as a = 5.791 A, b = 3.150 A measured 
with XRD. These lattice parameters are comparable to those 
reported for this boride in other disk alloys [23. 24], EDS and 
microprobe measurements on a sample aged 927 °C for 1,000 h 
(Figure. 5) indicated this phase had an approximate composition 
of 18.8Cr-15.8Mo-8.7W-4.5Ni-3.2Co-2.4Nb-2.4Ti-0.7Ta-3.8C- 
39. 7B in atomic percent. Table 2. Surveys in the literature have 
also reported M3B3 borides rich in Cr, Mo and W [1, 2, 23, 24], 
Mean equivalent radius and area fraction are tabulated versus 
sample condition in Table 2 for M3B2 borides. These borides 
were most numerous in the as-extruded sample, as fine particles 
less than 0.4 pm radius chiefly along the grain boundaries. The 
quantity observed along grain boundaries was greatly reduced by 
forging. Most borides dissolved after the supersolvus solution 
heat treatment, but widely scattered borides less than 0.5 pm 
radius reappeared both within grains and at grain boundaries for 
long aging times at 760 and 843 °C. The highest area fraction for 
each temperature was found at the longest aging time of 1,000 h. 
The overall highest area fraction and average size of these borides 
was observed for 1,000 h aging at 927 °C (Figure 5), where they 
grew into the largest minor phase particles found in this study, of 
1.19 pm mean radius. 

a Phase 

a phase particles had a tetragonal P42/mnm crystal structure, with 
lattice parameters a = 8.888 A, c = 4.593 A measured with XRD. 
Observed a precipitates were too small for composition 
measurements using the microprobe. Particles were separated by 
electrochemical phase extraction, and then analyzed by EDS in 
the SEM. This indicated a composition of 27.6Co-26.1Cr- 
14.2Mo-14.0W-l 1.7Ni-4.1Ta-2.4Ti in atomic percent, Table 3. 
The current measurements differ significantly from that reported 
in prior work, on superalloys having different overall 
compositions including lower Co content than LSHR [1, 2, 8, 9]. 

Mean equivalent radius and area fraction are shown versus sample 
condition in Table 2. o phase particles were found in a minority 
of examined cases, attaining a maximum area fraction of 0.00327 
at 843 °C for 1,000 h, (Figure 4). They appeared as irregular 
shapes at the grain boundaries, and seemed more elongated when 
nucleated in the grains. Some had distinct needle-like shapes. 
Several o phase particles were also found to be attached to small 
MC carbides. These a particles sometime appeared to have 
nucleated on the MC carbides, as this could require less energy. 
During aging at 843 and 927 °C, these a particles appeared to 
precipitate first at the triple points of grain boundaries, and later 
within the grains. However, they were not observed after 1,000 h 
at 927 °C, indicating they were not stable at this temperature. 
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Table 2. Measured minor phase mean equivalent radius (R, in nm) and fraction (f, in %). IG and GB represent intragrain and grain 
bou ndary locations, respectively. 


Sample ID 

MC j 

M 23 C 6 

1 m 3 b 2 ; 

1 7 1 

1 ® 1 

GB 

1 "! 1 

GB 

1 e 1 

GB 

All 

[ IG 

GB 

AN [ 

LbJ 

1 f 1 

R 


L_rJ 


R 

f 

LoJ 

1 f 1 

R | 

1 f 1 

R | 

m 

R 

f 

R 

f 

R 

f 

LAF-Forged 







- 

- 








- 

- 

- 

- 


LFF 





















LFW 

m 






□ 


m 











l_h 

LFW760-10 



m 


















LFW760-100 

m 















- 

- 

- 

- 


LFW760-1k 







ESQ 





BUfrll 



- 

- 

- 

- 

- 

- 

LFW 843- 10 







m 










m 


Wk 

■mail 

LFW 843- 100 









□ 





rzi 

163 

0.012 

216 

0.045 

196 

0.057| 

LFW843-1k 





















LFW927-10 

EH 














1 161 

0.006 

175 

0.053 

174 

0.059| 

LFW927-99 

EE23 


















IRI 


LFW927-1k 

IE3 






□ 








□ 







Table 3. Experi mentally measured versus thermodynamics 


predicted equilibrium minor phase compositions. 



Atomic % 

wm 

K3B 


WM 



MM 

mm 

MM 

c 

B 

wm 

LSHR Alloy Com 

positions 

wn 

EH 

■eh 

■a 


wm 

eh 

wn 

ED 

liffil 

■IBM 

gmi 


LFW760-1k 









wm 





StanDev 





eh 


EH 

EH 





Ni-Data 7 (1199 V) 




jHH 









Ni-Data 7 (760 V) 




MM 



eh 


EH 




PanNi 8 (1199V) 




■D 

wm 








PanNi 8 (978 V) 



MB 

1.3 

8.9 


EH 





wm 


LFW760-1k 




Kim 

EH 


■ EM 



EH 




StanDev 


■O 

EO 



EH 




ED 



Ni-Data 7 (760 V) 






eh 







PanNi 8 (760 V) 


mm 


MB 


ED 








LFW927-1k 




wm 

wm 

EH 

EO 

eh 

MM 

EH 

eh 

wm 


StanDev 


eh 

EH 

MB 

0.4 

0.1 

MM 

wm 

1.0 

1.3 

5.8 


Ni-Data 7 (927 V) 


•’■fiCT 







ED 




PanNi 8 (927 V) 


MB 







ED 





LFW843-1k 

Mean 



Km 

KH 



EH 

EH 

■ EM 




StanDev 


MB 

mn 

EH 


eh 

EH 

EH 

1.5 




Ni-Data 7 (843 V) 









EH 




PanNi 8 (843 V) 




EH 





■nm 





Modeling 

Equilibrium Modeling 

Using the CALPHAD multicomponent thermodynamics database 
Ni-Data 7 [25] and Thermo-Calc software [26], equilibrium phase 
fraction predictions versus temperature are shown in Figure 7. 
The Ni-Data 7 database was shown in our earlier work [18] to 
have better agreement with experiments in predicting y/y ' 
compositions and phase fractions than other databases. Note that 
this calculation was done using an effective alloy composition 
(second row of Table 3) obtained with a single equilibrium 
calculation at 1400 °C, before other secondary solid phase 
formation, to remove components related to primary particles, 
which are alumina M 2 0 3 (0.11 %) and M(C,N) (0.015 %). This 
alters the alloy aluminum and titanium concentration only 
slightly. The obtained effective compositions are used in all 
subsequent model calculations. Furthermore, the results in Figure 
7 were calculated without considering p and P phases, which were 
not observed experimentally. If they were considered in the 
calculation, p phase (with atomic composition of 
Coo. 30 Moo. 22 Cro. 17 Nio 15 Wo. 14 ) would become stable below 900 °C 
and reduce 0 phase stability. 


The MB 2 boride (basically TiB 2 ) is predicted to be stable between 
1070 to 1250 °C, as shown in Figure 7. However, the MB 2 phase 
was not found in any LSHR samples discussed earlier, probably 
due to a very low phase fraction and slow precipitation kinetics. 
The lower temperature M 3 B 7 Cr-Mo-rich borides are indeed found 
in the as-forged (LAF-Forged) samples (see Table 2), probably 
due to favorable precipitation kinetics over MB 2 particles. The 
subsequent 1199 °C solution treatment eliminated almost all of 
M,B 2 particles found in the as-forged condition. However, this 
was not compatible with a separate equilibrium calculation 
suspending the MB 2 phase, which suggested the stability of M,B 2 
extends up to 1245 °C. 

The MC and M 23 C 6 carbides that were observed experimentally 
are predicted by equilibrium calculations in Figure 7. MC 
carbides are predicted to exist at higher temperature, consistent 
with the experimental observation that MC exists in all LSHR 
samples. M 23 C 6 is predicted to be stable below 850 °C, which is 
qualitatively consistent with the experimental observations. 
Suspending a phase in the calculation increases the M 73 C 6 solvus 
temperature by 5 °C. Below 850 °C, formation of M 73 C 6 drives 
the dissolution of MC. This behavior is partially confirmed by 
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experimental observations summarized in Table 2, where total 
volume fraction of MC is reduced in 843 °C aging. This MC 
dissolution process is predicted to release titanium, which results 
in a slight increase of y ' phase fraction. 



Figure 7. Predictions of equilibrium phase fraction versus 
temperature, using Thermo-Calc with Ni-Data 7 database. 

Table 4. Predicted equilibrium phase fractions (contents) at the 
three aging temperatures. 


Tempera- 
ture, °C 

Phases 

Phase Fractions % 

Ni-Data 7 

PanNi 8 

760 

y’ 

53.6 

54.2 

MC 

0.048 

- 

M 23 C 6 

0.55 

0.65 

M 3 B 2 

0.38 

0.38 

Sigma 

6.8 

8.7 

843 

y’ 

50.7 

52.9 

MC 

0.17 

- 

M 23 C 6 

0.25 

0.61 

M 3 B 2 

0.38 

0.38 

Sigma 

1.3 

6.3 

927 

y’ 

46.4 

50.1 

MC 

0.28 

- 

M 23 C 6 

0 

0.56 

M 3 B 2 

0.38 

0.37 

Sigma 

- 

3.5 


The c phase is predicted to be stable below 860 °C and can reach 


a significant area fraction below 800 °C, according to equilibrium 
predictions. Omission of the boride phases in modeling increases 
a solvus temperature by 10 °C. The experimental observation of 
a phase at 843 °C aging confirms its stability, but at much lower 
observed area fractions due to slow precipitation kinetics. In 927 
°C aging, however, the initial observation of a particles cannot be 
explained with the equilibrium predictions, though final 
disappearance at 1,000 h aging does support the thermodynamic 
predictions. 

Table 4 summarizes the equilibrium phase fractions at the three 
aging temperatures predicted from both Ni-Data7/Thermo-Calc 
and PanNi8/Pandat, using the same assumptions discussed above. 
One major difference between the two databases is the transition 
temperature from MC to M 23 C 6 carbides, where Ni-Data 7 
estimates 850 °C and PanNi suggests a higher temperature of 984 
°C. Furthermore, PanNi 8 predicts a phase to be stable up to 
1016°C, much higher than that predicted by Ni-Data 7 database. 
However, from experimental data, it is not clear which database 
provides a better description of a phase stability. 

Along with measured minor phase compositions of 1 ,000 h aging 
samples, the predicted compositions are also listed in Table 3 
using both Ni-Data7/Thermo-Calc and PanNi8/Pandat at the 
corresponding aging temperatures. For MC carbides, since they 
pre-exist before the aging treatment, the predicted equilibrium 
compositions at the 1199 °C solution treatment are also shown in 
the table. The measured MC compositions in sample LFW760-lk 
can be largely explained by starting with the Ni-Data 7 predicted 
1199 °C composition as a baseline, with small modifications 
towards that predicted at 760 °C. Since MC is not stable at 760 
°C by the PanNi 8 database, the MC composition is listed at its 
lowest stable temperature of 978 °C. The PanNi 8 composition 
underestimates the Nb and Ta contents, while overestimating Cr 
and Ti contents. 

For M23C6 carbides, there are large discrepancies of the predicted 
compositions from both databases versus the measured 
compositions. PanNi 8 predicts nearly pure Cr 23 C 6 carbide 
composition, with only minor Co substitution for Cr. Ni-Data 7 
predicts additional Mo and Ni substitution for Cr. Experimentally 
measured compositions suggest much higher Co, Ni and W 
contents, and a significantly lower Cr content. 



Figure 8. Trajectory evaluation of carbides at the three isothermal aging temperatures. Driving force of MC (M 23 C 6 ) is shown in the X- 
(Y-) axis. The stars represent the experimentally observed carbide fractions. The precipitation trajectory is formed by connecting the stars. 
The solid and dashed curves represent constant phase fractions curves of MC and M 23 C 6 carbides, respectively. 
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Predictions of M 3 B 2 compositions from both databases are very 
similar. However, both predictions overestimate Mo content, and 
underestimate W, Ni, Nb and Ti. The predicted 0 phase 
compositions underpredict amounts of Mo, Ta and W, and 
overpredict Cr content. PanNi 8 does predict a closer W content 
than does Ni-Data 7. 

The discrepancies between model predictions and experimental 
observations could be caused by either the inaccuracy in 
thermodynamic databases or the observed phases not representing 
equilibrium conditions. If thermodynamic databases are indeed 
accurate, three explanations concerning non-equilibrium 
partitioning are listed below: 

• Partial solute trapping and para-equilibrium constraint 
might explain the underestimation of particle compositions 
for slow diffusing elements. For example, the observed 
lower W partitioning ratio between y ' and y suggested 
partial solute trapping [18]. For minor phases studied 
here, observed lower Mo contents in M 23 C6 and M 3 B 2 
might be driven by partial partitioning. However, the 
opposite effect is observed here for W in M 23 C 6 , M 3 B 2 and 
0 particles, where actual W contents are higher than the 
equilibrium predictions 1 from CALPHAD thermodynamic 
databases. This could suggest a need to refine the 
multicomponent databases. 

• Formation of M 23 C6, M 3 B 2 and 0 particles mainly occurs 
on grain boundaries. The deviation of particle 
compositions can be driven by the enhanced diffusion on 
grain boundaries and multicomponent flux balance 
requirement. 

• Compositional segregation to the grain boundaries was 
observed in Ni-base alloys. For example, substitutional 
elements, Mo, Nb and W were shown to partition to grain 
boundary [27, 28], This could lead to very different local 
compositions for the formation of grain boundary 
particles. This might also explain the initial precipitation 
of 0 and M 23 C 6 at grain boundaries for 927 °C aging. 

Carbide Precipitation Trajectory Analysis 

The three plots in Figure 8 summarize the precipitation trajectory 
of MC and M 23 C 6 carbides based on an unstable equilibrium 
fraction- supersaturation relationship where the precipitation 
driving force is balanced by capillary energy. This calculation 
was performed using Thermo-Calc with the Ni-Data 7 database, 
assuming y' at equilibrium, and no M 3 B 2 and 0 phases present. 
This assumption is reasonable, as y r precipitation kinetics is faster 
than for the carbides, while that for the boride and 0 phases 
appears to be even slower. The solid curves for MC and dashed 
curves for M 23 C 6 carbides represent constant phase fractions 
within the X- (Y-) axis of driving forces of MC (M 23 Cs). At the 
origins of the plot, driving forces for both carbides are zero, 
representing final equilibrium conditions. The star data points 
represent the experimentally observed total carbide fractions at 
different conditions. A smooth path would ideally be expected 
from the starting as-quenched (LFW) condition to the final 1,000 
h aging condition, moving toward the origin. However, for both 
843 and 927 °C aging treatments, the 100 h conditions seem to 
move away from equilibrium. Furthermore, due to slower 
kinetics, the condition of 760 °C for 1,000 h is still not close to 


1 Predicted M 23 C 6 , M 3 B 2 and 0 compositions under supersaturated 
unstable equilibrium condition still have the same discrepancy versus 
experimental measurement. 


the equilibrium point. This information helps us understand the 
precipitation path along the thermodynamic supersaturation, and 
provides calibration input for the multiphase simulations 
discussed next. 

Precipitation Modeling for Carbides 

PrecipiCalc software has been successfully applied to capture the 
location-specific y' microstructure’s evolution in a y matrix for 
Ni-base disk alloys during non-isothermal processing [17, 18]. In 
this work, PrecipiCalc software was used to simulate the 
concurrent y ', MC and M 23 C6 carbide microstructural evolution 
during isothermal aging, in order to assess the feasibility of 
modeling complex precipitation interactions and to identify 
potential model limits. Since PrecipiCalc utilizes CALPHAD 
thermodynamics and mobility databases, we can then also assess 
the impact of the discrepancies in phase fraction and compositions 
discussed earlier on the precipitation predictions. 

In order to achieve a manageable simulation, the focus is placed 
on the carbides, assuming both boride and 0 phases can be 
suspended. This is reasonable as their formation kinetics is 
slower than carbides, as shown in Table 2. Even though M 3 B 2 
achieves 0.8% phase content at grain boundaries in aging 927 °C 
for 1,000 h, its impact on carbide compositions is still judged 
small. The y' phase reaches equilibrium at a much faster rate than 
the carbides. 



Figure 9. Predicted secondary y' mean size evolution in LSHR 
from PrecipiCalc (curves) [29]. The predicted size in as-quench 
(LFW) sample is shown at 0 h aging time. 

According to PrecipiCalc simulations, within one hour of aging 
above 760 °C, y' volume fraction should be at equilibrium. To 
further simplify the simulation, one could also consider removing 
the equilibrium y' phase fraction and compositions to focus on the 
carbides’ precipitation in the y matrix. However, in this 
isothermal aging study the y ' phase was included, using the 
parameters found in the previous study [ 1 8] including a starting 
microstructure of 14 nm mean radius and 47 % phase content 
predicted from a separate PrecipiCalc simulation following the 
fast water quench (LFW) process. Results of predicted y ' size 
evolution at the three aging temperatures are shown in Figure 9, 
and further details can be found in [29], After aging the grain 
boundary MC carbides constituted no more than a quarter of the 
total MC phase fraction, according to Table 2. Hence, MC 
carbides are modeled as intragrain particles with a homogeneous 
nucleation model, while later comparisons are made to total 
experimental MC size and fraction measurements. For M 23 C6 
carbides, PrecipiCalc' s classical heterogeneous grain boundary 
nucleation model for incoherent particles is used, with a 25 pm y 
grain size and a wetting angle of 33°, roughly representing the 
mid-value of the observed M 23 Cg morphology, which did often 
vary widely. PrecipiCalc' s 3-dimensional volume diffusion- 
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based growth model is utilized for both MC and M23C6 carbides. 
Applicability of this growth model to grain boundary M23C6 
carbides will be discussed later. There is a prototype of 
PrecipiCalc with a revised grain boundary collector plate growth 
mechanism [31] which could also be relevant for M23C6, but the 
current implementation is limited to only one precipitate phase, so 
it is not used here. 

The PrecipiCalc simulations discussed below rely on pre- 
processed PrecipiCalc parameters based on Ni-Data 7 and the 
NIST mobility database [30] using the unstable equilibrium 
approximation of the previous trajectory model. The pre- 
processed parameters help computational efficiency and also 
improve simulation stability. Once these parameters are 
generated, subsequent PrecipiCalc simulations no longer require 
CALPHAD free energy optimizations, which can be challenging 
for multi-phase precipitation simulation. 

The starting microstructure at the beginning of the isothermal 
PrecipiCalc simulations is based on the experimental results for 
the LFW forged and water quenched sample. That is, (1) no 
M23C6 particles are assumed present; and (2) the MC particles are 
assumed to be normally distributed, have a mean radius 2 of 173 
nm and a standard deviation of 60 nm. Total MC phase content is 
taken as 0.164 % 3 . 

In order to model the observed experimental data, calibration of a 
set of precipitation kinetics parameters was first necessary. The 
objective here was not to just model experimental results as 
accurately as possible, but to search for a minimum set of physical 
parameters to capture the observed behavior, in order to assess the 
sources of model inaccuracy. This calibration determined the 
following parameters. 

• Carbide/y interfacial energy was found to be 0.8 J/m 2 for 
both MC and M23C6 at all three temperatures. This higher 
than expected interfacial energy for intragrain MC 
carbides is probably necessary to account for the under- 
estimation of slow diffusing Ta in MC (see Table 3). 

• A diffusion scaling factor of 0. 1 for M23C6 grain boundary 
carbides was needed and is justified as following. Though 
grain boundaries offer a faster diffusion path, volume 
diffusion of solute through the grains’ microstructure 
containing a high fraction of y' particles could dominate 
the formation rate. The reduction of diffusivity for M23C6 
carbides is needed to correct the y' diffusion barrier. In 
addition, the thermodynamic model underestimates the 
amount of the slow diffusing element W in M23C6 (Table 
3), and results in an overestimated particle growth rate. 
Scaling down diffusion is thus needed to move toward the 
correct growth rate. This can also be related to an 


2 Small inaccuracy introduced here due to treating experimentally 
observed 2D particle sizes as 3D sizes to compare with 
PrecipiCalc simulations. 

1 The MC carbides in the as-quench sample are responsible for 
pinning the y matrix grain boundaries. With the measured MC 
size and phase content, the Zener pinning model, based on 
spatially uniform distribution of pinning particles, predicts a y 
grain size of 16 pm, under-estimating the measured grain size of 
25 pm. This might suggest a spatially non-uniform MC carbide 
distribution. No abnormal grain growth was observed after the 
supersolvus heat treatment applied here. 


incorrect use of a 3-dimensional growth model for grain 
boundary precipitation discussed later. 

The results of the simulation are summarized and compared with 
experimental data in several plots within Figure 10. Time 
evolution of experimental mean radii, phase contents, and derived 
number densities are shown as discrete points based on 
experimental data in Table 2. PrecipiCalc results (labeled with 
“PpC”) are shown as thick curves for mean size and fraction 
evolution. The microstructure of as-quench (LFW) sample is 
shown at the 0 h aging time. 

Driven by the underlying thermodynamics, MC dissolution and 
M23C6 growth are predicted at 760 and 843 °C, though there is no 
clear experimental evidence of MC dissolution at 760 °C. Overall 
qualitative trends of carbide precipitation are captured by 
PrecipiCalc, but there are significant quantitative discrepancies. 
In addition to the incorrect carbide compositions from bulk 
thermodynamics discussed earlier, a large part of the 
discrepancies can be rationalized by the use of (1) the particle 
growth model based on 3 -dimensional volume diffusion, (2) the 
single set of nucleation parameters on grain boundary, and (3 ) the 
mean field approximation where the carbides are distributed non- 
uniformly. The mean field assumption of PrecipiCalc considers 
all particles as surrounded with a mean far-field composition. 
Individual particles can interact with any other particles 
instantaneously. In reality, the carbides in LSHR depart from the 
mean field approximation. For example, grain boundary M13C6 
carbides interact with only a small portion of MC carbides on or 
near the grain boundary. This also means the local surrounding 
compositions for the particles can strongly depend on the 
locations of carbides. As a result, use of the mean field 
assumption in this PrecipiCalc simulation could overestimate the 
communication between the carbides. This can explain the faster 
kinetics of dissolution and growth of carbides predicted by 
PrecipiCalc. 

An important parameter of precipitate evolution during growth or 
coarsening is the time exponent for mean particle sizes. For 3- 
dimensional diffusional growth or coarsening of intragrain 
particles such as y’, the time exponent is 1/2 or 1/3. For grain 
boundary particles having no interactions with intragrain particles, 
experimental evidence and modeling have shown a lower 
coarsening time exponent, down to as low as 1/5 [32], If one 
includes interactions with intragrain particles, growth and 
coarsening of grain boundary particles can be slowed down to a 
time exponent of 1/5 and 1/6 [31], respectively. As shown in 
Figure 10. the evolution of M23Q mean size at 760 and 843 °C, 
follows an exponent of 1/3.5 and 1/3.7, while PrecipiCalc follows 
a higher exponent. From the derived number density plots, it is 
shown that the experimental number density is increasing at 760 
°C and following a time exponent of 1/3.5, suggesting the 
nucleation process is ongoing up to 1,000 h. At 843 °C, 
nucleation of M23C6 seems to stop at 100 h. Hence, actual M23C6 
precipitation includes both nucleation and growth mechanisms. 
PrecipiCalc’ s nucleation kinetics is currently too fast and sharp, 
likely caused by the use of a single set of nucleation parameters, 
such as the wetting angle. In reality, a distributed and smooth 
nucleation behavior due to variation of grain boundary 
misorientation might be needed. 

The mean sizes of experimental MC carbides are found to be very 
stable. By analyzing the data in Table 2, the time exponents of 
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intragrain and grain boundary MC carbide mean radii are within 
1/20 to 1/60. The slow kinetics of MC carbides is likely caused by 
their spatially non-uniform interaction with the M 23 C6 carbides 
and diffusion barriers related to the large volume fraction of y' 
particles. 

Most existing precipitation models focus on modeling single 
phase grain boundary or intragrain particles. The carbides in 
LSHR show complex precipitation behavior within a dual phase y 
/ y ' matrix, and suggest trends beyond current precipitation 
modeling capabilities. Note that these modeling challenges and 
limitations discussed above for carbides also apply to the o phase 
precipitation, Further model extensions appear necessary to fully 
describe the carbide precipitation behavior observed in LSHR Ni- 
base superalloys. After refinement of these carbide precipitation 
simulations, the lessons learned can then be applied to include 
simultaneous precipitation of the a phase. 

Summary and Conclusions 

The predominant minor phases in disk superalloy LSHR were 
characterized after varied heat treatments. The results were used 
to assess existing model capabilities for predicting equilibrium 
states and precipitation for these multiple phases. It can be 
concluded from this work that: 


1) Initial trials indicate joint precipitation of y', MC carbide, 
and M 23 C 6 carbides can be simulated using models such as 
PrecipiCalc. 

2) Current model limitations are identified and additional model 
extensions have good potential to accurately describe 
precipitation of y ' and carbides during aging heat treatments. 

3) This may require modifications and extensions to current 
treatments of fundamental multicomponent thermodynamics 
and multi-phase precipitation models. In particular, the focus 
should be placed on elemental partitioning, grain boundary 
diffusion, segregation and distributive nucleation, and the 
interaction of spatially non-uniform particles. 
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Figure 10. Comparison of LSHR carbide microstructure evolution between experiments (points) and PrecipiCalc predictions (thick 
curves labeled with “PpC”). Both MC (squares and solid curves) and M 2 3 C 6 (diamonds and dash curves) are shown. The as-quenched 
carbide microstructure is shown at 0 h with hollow symbols. Thin lines represent power law fits to the experimental data, suggesting the 
time evolution exponent. The number density is calculated as phase fraction divided by particle volume at mean size. 
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